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TECHNICAL NOTE NO. 1013 

VELOCITY DISTRIBUTION ON WING SECTIONS OJ ARBITRARY SHAPE 
IN COMPRESSIBLE POTENTIAL PLOW 
II - SUBSONIC SYMMETRIC ADIABATIC PLOWS 
By Lipman Bers 

SUMMARY 


This paper extends the method of computing the pressure 
distribution along a symmertioal profile of arbitrary shape 
given in NACA Technical Note No. 1006 under the assumption 
of a linearized pressure-volume relation to the case of an 
everywhere subsonic flow satisfying the rigorous adiabatic 
equation of state. Either the stream Mach number or the 
maximum local Mach number may be prescribed. 

The actual applicability of the method depends upon the 
development of efficient procedures for numerical integra- 
tion of linear partial differential equations with variable 
coefficients , 

Tables and graphs of functions needed for the computa- 
tion are given, and a velocity correction formula proposed 
by Garrick and Kaplan is discussed. 


INTRODUCTION 


In a recent paper (reference l) the author presented an 
effective method of computing the velocity distribution along 
a symmetrical profile of arbitrary given shape under the as- 
sumption that the potential compressible flow past the pro- 
file obeys the so-called linearized pressure-density relation. 
In this paper this method is extended to the case of the ad- 
iabatic preB sure-density relation which actual flows satisfy 
with a high degree of accuracy. In fact, the method presented 



here ia applicable to any empirically or theoretically given 
pressure-density relation, provided the flow remains every- 
where subsonic. 

The computational labor involved in applying this method 
is considerably greater than in the case of the linearized 
pressure-density relation. The aotual computations could 
hardly be attempted without using either automatic computing 
maohines or electrical models. In fact, at the present state 
of mathematical analysis, it is doubtful whether reliable 
theoretical results on pressure distributions for high sub- 
sonic speeds and on the values of the critical Mach number 
for various profiles could be obtained by any method without 
extensive numerical computations. On the other hand, the 
present ideas concerning what constitutes difficult computa- 
tions will doubtless undergo a very radical change in the near 
future due to the development of modern automatic computing 
machines . 

In the author 1 a opinion, the merit of the present method 
lies in the fact that it eliminates the mathematical diffi- 
culties peculiar to the fluid dynamical problem considered 
and reduces the solution of thiB problem to the numerical 
integration of a linear partial differential equation, This 
latter problem is of prime importance to all branches of ap- 
plied mathematics and is the subject of many investigations. 

The present method is based upon the application of a 
transformation which for compressible flows takes the place 
of conformal mapping. This transformation was introduced in 
a previous paper. TSee reference 3.) The pressure distribu- 
tion problem is first reduced to the determination of a func- 
tion f(uj) defining a point-to-point correspondence between 
the prescribed profile and a circle. Then it is shown that 
this function satisfies a functional equation of the form 
f = F(f) where F is a certain functional operator. In 
order to determine f, a trial solution f 0 is chosen, and 
the functions f a = f(f 0 ), f 3 » fC^),.., are computed. On 

the basis of computations carried out for a simple special 
case (V = -l), it may be expected that the sequence f lT 
fa,... rapidly converges toward the desired function f (uj) . 
The computation of the transform F(f n ) involves an integra- 
tion of a linear partial differential equation previously 
mentioned, 

On the basis of the theoretical results obtained, a dis- 
cussion of a "velocity correction formula" proposed by 
Garrick and Kaplan is given. 
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SYMBOLS 


constants (appendix I) 

auxiliary functions defined by equation (36) 
speed of sound 
positive constants 
non-Buclidean length element 
E(P) domain exterior to the profile P 

J[f(u>« ),u>] functional transformation defined in section 8 
f (uj) 


A, B, b 

cons 

a(£), b U) 

auxi 

a 

spee 

o 

o 

posi 

dS 

non— 


Sik 

HU) 

h(uj) 

M 


function defining the mapping of the circle into 
the profile 

coefficients of the metric generated by the flow 

harmonic function defined by equations (42) and 
(43) 

function defined by equation (40) 

Mach number 


P 


P 


profile in the z-plane 

lattice points in the Z-plane (appendix I) 
pressure 
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ft 

fti 

q. 

< 1 * 

r 

S 

a 

S 


u, v 

U* , V* 

w 

V* 

*, y 
x, y 

z 

z i» z T 
Z 


a 



c 

0 


dimensionless 

^peed 

<l/a 0 

as funotion 

of q * a 

dimensi onless 

speed 

q/ a o 

as function 

of q* 


speed 

distorted speed 

modulus of £ 

length of the curve P 

arc length measured along the profile 
function defined by equation (4) 
velocity components 

components of the distorted velocity 
complex velocity 
distorted velocity 

Cartesian coordinates in the z~plane 

Cartesian coordinates in the Z-plane 

complex variable in the plane of the flow 

leading and trailing edges 

auxiliary ooraplex variable 

angle at the trailing edge 

coefficients of the difference equation 
(appendix I) 

coefficients of the Liebmann transformation 
exponent in the pressure-density relation 
complex variable in the plane of the circle 
slope of the profile P 

angle between the velocity vector and the x-axis 


0 



a t> 
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MmJ 

p 

a 

T 

X 

cp 

$ 

* 

uj 

( ) 

( ) 

( ) 

( ) 
C~) 


value of < 1 * for M = M B 

Cartesian coordinates in the £-plane 

density 

dimensionless length parameter or the profile 
argument of £ 

auxiliary function defined hy equation (35) 
velocity potential 

function proportional to the velocity potential 
function proportional to the stream function 
argument of a point on the circle | £ | » 1 


X 

partial derivative 
(similarly for 

S o 
- Hi 

( ) with 

£, etc. 

0 

value of ( ) 

at 

a 

stagnation p 

CO 

value of ( ) 

at 

infinity 

max 

maximum of ( 

) 

i n 

the domain c 


boundary value 

of 

( 

) 


ANALYSIS 

1. Basic Relations 


A steady two-dimensional potential flow of a compressi- 
ble fluid is characterized by the quantities: p (pressure), 

(density), q (local speed), and 0 (angle between the 
elocity vector and some fixed direction, say the direction 
of the x— axis). The components of the velocity vector are 
given by 

u = q cos 8, v = q sin 8 


the complex quantity 
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W as U - lV 


4« 


-16 


is called the (conjugate) complex velocity. Pressure and 
density are connected by the adiabatic relation 



P 0 /P 


Y 

o 


where the subscript 
point. The exponent 
Instead of the speed 
sionless quantities 


o refers to the fluid at a stagnation 
Y is a constant; for air Y = 1.405, 
q it is convenient to use the dimen- 


and M 
a o 


q. 

a 


where 



is the speed of sound. 


Bernoulli's theorem implies the relations 


i 



( 1 ) 


( 2 ) 


It will be convenient to use the following two func- 
tions q* and T of the speed which are defined in the 
subsonic range by 
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and 


Here 



(3) 

(4) 


is the critical speed for which M = 1, A simple computation 
yields 



f y-1 q a 7 Y+l q a 

J_ ^ a O^ j ^ a Q 2 

Y -i q. a , 1 y+i g a 

2 a 0 a v 2 a 0 a 



T 

v V 1 ~ 


-^1 


Y-l q’ 


1 - 


Y+l q' 


iPFS-/- 


(5) 


Y+l q" 




J 


where 



( 6 ) 


It is clear that the "distorted speed" q * is an increasing 
function of the dimensionless quantity q/a Q and that 


q* = 0 for q = 0, q* = 1 for q = q s 
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q = a Q Q(q* 2 ) T • T(q* 2 ) • (7) 

The functions discussed In this section are tabulated 
(for X = l.k05) in tables 1 to k- and plotted in figures 1 
to 3. Tables 1 and 2 on which the other computations were 
based were taken over free, tables computed on the I. B. M. 
Automatic Sequence Controlled Calculator. (See reference 3*) 

The complex function 


w* = q* e 


-IS 


( 8 ) 


is called the distorted velocity. The knowledge of w* im- 
plies that of w, since 




Remark 1 : For an isothermal flow (X = 1) the preced- 

ing formulas must be replaced by the following 


P » P Q e 


-q 2 / 2a Q 2 -q 2 / 2a 2 

0 , P « Po© 0 


M = q/a 


t* 


j/fa 


J 1 - * 2 / a 0 2 




v 


T = /i - *** 


a o 2 


Rem ark 2: The preceding formulas, with the exception of 
those for p, remain valid if the pressure -density relation 
I 3 assumed in the slightly more general form 
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t 


p = A + B 

where A and B are constants. 

Remark S t The method described in this report can he 
applied to any pressure-density relation, provided the 
flow remains subsonic. 


2. The Boundary Value Problem 

Consider a symmetrical profile P in the plane of the 
complex variable z = x + iy . It will be assumed that the 
x-axis coincides with the axis of symmetry of the profile and 
that the trailing edgo zj is to the right. Let S be the 

length of the curve P and s the arc length measured along 
this ourve from the point zj in the counterclockwise direc- 
tion. The points on the profile may be characterized by the 
dimensionless parameter 




* 


cr ■ 2 tt -| (9) 

Thus ct = 0 corresponds to Zrp and a = tt to the leading 
edge z . Furthermore, let ©(a) denote the slope of P 

Ju 

at a point corresponding to ct, This angle is defined aB 
the positive angle by which the positive x-axis must be 
turned to move it parallel to the tangent vector to P point- 
ing in the direction of increasing ct. (See fig. 6.) Thus 
the function ©(ct) depends only upon the shape of P and 


0(0) s tt - a/2, ©(tt) = 3tt/2, ®(2tt) = 2tt + a/2 (lO) 

where a is the angle at the trailing edge ( 0 ^ a ^ tt - ). 
Furthermore, 


q( 2tt - ct) = 3tt - 0(a) (ll) 

The equation of the profile P can be written in the form 




z = 


Zip + S 


o 


0 - a 5 2tt 


(13) 
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Consider now a steady subsonic circulation-free poten- 
tial flow of a compressible fluid past P. Let a 0 cp(x,y) 

"be the velocity potential. Then 


u = a 0 


dcp , dcp 

— and v = ar — 


&x 6y 

and the magnitude of the velocity vector q is given "by 


(13) 


q. = a 0 

so that the continuity equation 


Mf : W 


(14) 


1_ (£- 


^p Q bxj Sy \j 


!V o 

dy VPo 


(15) 


■becomes a nonlinear partial differential equation for cp. 
In fact, this equation may "be written in the form 




Y - 1 


V) 


29 cp cp 
^x y xy 


+ o- 


Y - 1 a 

9 X 


Y + 1 




') 


^yy " 0 


(16) 


where subscripts denote partial differentiation. Por the 
special case Y = -1, this becomes the equation of a minimal 
surface. 

The velocity potential cp(x,y) is determined as the 
(one-valued) solution of this differential equation satisfy- 
ing certain auxiliary conditions. The first of these states 
that the undisturbed flow has the positive x-direction. It 
may be written in the form 


lim 


OD 


Sep 

sf > °' 


lim = 0 

s y 


(17) 
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Since 




(18) 


u 


this condition is equivalent to the following 

lim 0=0 

z CO 


(19) 


The second condition is a "boundary condition expressing 
the fact that the profile P is a streamline! 


dn 


on 


( 20 ) 


where d/dn denotes differentiation in the normal direction, 
In view of the symmetry of the flow this condition may "be 
restated as follows 


0 = 



on the 



"bank of 


( 21 ) 


Finally, a condition is needed which determines the 
magnitude of the speeds involved. This can "be done in two 


ways. It is possible to prescribe the speed 
undisturbed flow. The condition then reads 


of the 


or 



% 00 
a o 


( 22 ) 


lim q = q (23) 

CO 

Z CO 


Alternatively, it is possible to prescribe the maximum local 
Bpeeds q ffiaX . The condition reads 


max 



(24) 
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or 


max q 


^■max 


(25) 


In the first case, the stream Mach number is prescribed, 

in the second, the maximum local Mach number M max , In the 
second case, the condition imposed on cp has the character 
of a boundary condition, since the maximum speed is attained 
at the boundary. (The proof of this statement will be found 
in reference 4; it also follows from the considerations of 
s ec . 4. ) 

In the following, this boundary value problem will be 
reduced to a mapping problem. 


3, Mapping of the Profile 


into a Circle 


The hodograph of 
domain in the w-plane 
35(P) exterior to P 


the compressible flow around P is the 
(w = u - iv) into which the domain 
is taken by the transformation 


w = 




(26) 


In the case of a symmetrical flow the hodograph is a doubly 
covered Riemann surface with a branch point at w = q^ as 
shown in figure 5, This surface is transformed into the "dis 
torted hodograph" by the transformation 


w* = q*(|wj/a 0 ) (27) 

where q* is the function of q = | w J defined by (3). 
Finally, the distorted hodograph is mapped conformally into 
the domain I £! > 1 of the auxiliary £-plane by an analytic 
function ^ C w ^ ) satisfying the conditions 


i = 


a t 

dw* 


> 0 for 


w * = q * 


Co 


= ^(^/a 0 ) 


In this way a transformation 


I = ! (x,y ) , ^ = *l(x,y ) 


(28) 
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of the domain E(P) into the domain |£| > 1 is obtained. 
It is seen that the points 


Z = °= , Z = Z ^ , Z = Z rp 

are taken into the points 

l = ®, £ - -l, £ - l 

respectively, and that the horizontal direction at infinity 
is preserved by the mapping! 


11a > 0, 

2 — >* 03 


lim 

2 03 


at 

dy 


0, lim ~- 
a — os ox 


0, 


lim 

a -*■ co 



0 


(29) 


The resulting mapping of the profile P into the unit 
oirole can be described by means of a function 

a m { (cu ) (30) 


such that a point z of P corresponding to the parameter 


value c = f (in) is taken into the point £ ■ e 
tion (30) is an increasing function and 

f(0) = 0, f ( 2tt ) = 2tt 

furthermore, by reasons of symmetry 

f(2u - a) = Stt - f(a) 

so that 

f (rr ) * tt 


iin 


The func- 


(31) 


(32) 

(33) 


It will be shown that the knowledge of the function f(cu) 
implies that of the velocity and pressure distribution along 
P. 

Remark 1 1 It has been shown (reference 2) that the 
transformation (28) is conformal with respect to the following 
Riemann metric! 
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dS® 


gn dx® + 2g ia <ixdy + gsa<iy a 


(34) 


where 


g 13 . = 1 - M 8 sin 8 6 

3 

g ia = M ain 8 sob 8 > 

g aa = 1 - M 8 cos 3 6 


(35) 


It is instructive to consider three special oases. Tor 
an infinitely slow flow (i.e 4 , for an inc ompr es s ihle fluid) 

M = 0 and dS a = dx 3 + dy a . Then the mapping just constructed 
is the standard conformal mapping of a profile into a circle, 
lor a uniform flow in the x-direction 0 = 0, M is constant 
and dS 3 * dx® + (l - M a )dy®. The introduction of the metric 
(34) is equivalent to the well known Prandtl-G-lauer t contrao* 
tion. Finally, in the case Y * -1, the mapping (28) coin- 
cides with the mapping defined hy the equation 


« 


constant 


0 + t) = ^ + iV{/ 


where ^ is the stream function of the flow, (See reference 

1 . ) 


4, Computation of the Velocity in the £-Plane 

From the way in which the mapping of the z-plane into 
the J-plane has "been defined it follows that the distorted 
velocity 

* * -ie 

w* * q*e 

is a one-valued regular analytic function of 5 for |t| > 1. 

It vanishes only for t '=* -1 and (unless a = 0) for £ ** 1. 
Moreover, the funotion 


w *£ i + a/tr 

< £ + lX £ - l) a/TI 


x(£) = 
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is continuous for ! £ | ^ 1 and is nowhere in this domain 
equal to either zero or infinity. (For the proof of this 
last assertion } see reference 1, sec. 3.) 


Set 


log X B A + iB 


(36) 


Then A + iB and hence B - iA are regular one-valued ana- 
lytic functions. Hence A(£) may he computed in terms of 
the boundary values of the function B(£). For j £| > 1 
this is done by means of the P oi s s on-Schwarz integral formula 


B(£) 


iA(U 



i<Ju. , 
e )du> 


iA( c») 


which implies that 


a(£) 


W r sln(T - Ul) A („), 

n J l-2r cos(T-Uj)+r 3 

o 


( £ * re iT ) 


For |£| « 1 the formula to be used reads 

A(e luj ) = - ^ ^ ^B(e ilW ' iT ) - B(e iu> “ 1T )j> cot 1 dT + A(a») 

0 

(See, for instance, reference 5, p. 243.) Of course, it is 
also possible to determine first A(e^ a) ) and then to compute 
A(re*‘), r > 1 by the formula 


A(re* T ) 



1 - r a 

2r cos(t - uj) + r® 


A(e iU) )dU> 


If the value of 8 at a point of _P corresponding to 
the parameter value a is denoted by 8(a), 
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B(e iu 0 

and 


- SCf (u»)3 - arg (l + e icu ) - & arg(e iU) - l) + ^1 + (S?) 


erg (l ♦ .»> 

U)/ 2 + TT 


for 0 < UJ < TT 

for tt < u> < 2tt 


arg(e iuJ - 1) - HI + 1 
2 2 

^furthermore , by virtue of (36) and (21 ) 


^ f®[f(u>)] - n for 0 < to < tt 

e [f(uO] * < y 

^0[f (hJ) ] - 2 tt for TT < U) < 2 tt J 


Also 


log q* * A + log 


(£ » D(j - 1 ) 

— a iMieiniUi a. i m ~ — 

£ i+a/TT 


a 

IT 


( 38 ) 


and if the value of q* at a point of P belonging to the 
parameter- value a is denoted by q*(a) , then 


log- (to) ] ■= A(e iW ) + 1 


og 


r »+& 

<12 17 


C 0 3 


CU 


sin 


*> 


By noting that A(«) equals the logarithm of the value 
of the distorted speed at infinity and combining the preced- 
ing formulas , the following expressions for q* are obtained. 


?* £<“>)] = 5 * 2 




where 


cos Hi 


sin — 

2 | 


2 


a . 

117 h (0) ) 

© 


(39) 


.TT 


^ J j©[f(u)+ t) 3 - §[f(.UJ" t)]J>. oot J 4T (40) 
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and 


<£) 


<t + D(t- i) 




t 


l+a/rr 


ML) 


(41) 


where the (harmonic) function H(£) is given by either of the 
two equivalent formulas: 


H 


/ an 

r »ln(r - m) ( a[f . jga. uAiu, ( 42 ) 

1 - Sr cos(T-iu)+r s L 3n J 


or 


aTT 


H(£) « - i- 


1 - r' 


2tt l~2r o 08 (t-tu)+r a 
o 


h(uj)du) (43) 


In both cases £ = re^ T . 


If in the original boundary value problem the value of 
the speed at infinity has been prescribed, then q*^ is a 

known constant. If, however, the value of the maximum local 
speed is prescribed, then q* ie determined from the rela- 
tion ® 


2 




sin 




(44) 


where q* denotes the (given) maximum local distorted speed, 
max 

The funotion 8 also is uniquely determined by f(cu), 
for being a harmonic funotion of | and n and satisfying 
conditions (38) it is given by the formula 


e(C) 



1 - r a 

Sr cos ( T - uj) + r s 


{■ 


0[f (U))] 



dtu 



an 


1 - r“ 


l~2r oos(T-u))+r* 


|e[f(ch)] - 2 ttJ, dcu (45) 



NACA TN No. 1013 


18 


It should "be observed that 

lim 6=0 (46) 

£ — > <* 

No difficulties are involved in computing the above- 
mentioned functions since the integral defining the function 
h(in) is a proper Niemann integral* In fact, for T = 0 the 
integrand is equal to 

40* [f (u>)3 f* (u)> - iiSLtZl 

TT 

In computing this integral the expression in th e braces must 
be treated as a periodic function of tu with the period 2tt. 

After ci 181 and 6 have been determined the quantities 
q/a 0 , T, p, p, and gj.^ can be computed as functions of 

i and r\ . In particular, the boundary values q(cr) of q 
are given by 

q [f(u))3 = a 0 Q, <{V [f(o))] 3 | 

Thus the statement that the knowledge of the function cr = f(uu) 
yields the pressure distribution along N- is verified. 

In the preceding considerations no properties of the 
function f(u)) except the boundary conditions (31) were used. 
Therefore the following statement is valid. 

To every function f(uu) (which is not necessarily the 
mapping function) satisfying the conditions f(0) = 0 and 
f(2Tr) = 2tt there corresponds an analytio function w* = q*e“^° 
defined by equations (39) to (45). This function is uniquely 
determined by f (tu) , the function 0 (ct) which depends only 
upon the shape of the profile P and either of the two 
parameters q* and q* which are connected by relation 
(44). max 

If jw*| < 1 everywhere, then the functions q/a 0 and 

T also can be computed as functions of £. This will always 
be the case if the prescribed value of q* ax < ~ • 

If the function f(uu) also satisfies condition (31) (as 
the actual mapping function does), then the analytic function 
w* will satisfy the condition 
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V* (£) = v * tp (4?) 

(The "bar denotes the conjugate complex quantity.) Ia this 
case the function T = T(|,n) corresponding to f(u>) will 
possess a symmetry with respect to the |~axisi 


E( t.-n ) = T( |,ri) 


(48) 


In the following condition (3l) will always he assumed. 


5, Computation of the Potential in the £-Plane 

It is known mat the potential cp considered as a func- 
tion of u* and v* (w* = u * - iv*) satisfies the linear 
partial differential equation 

-A- (k (L IseA = o (49) 

du* \T Su V dv* \T dv*/ 

(See reference 2, p, lip.; This elliptic equation is invariant 
with respect to conformal transformation so that <p considered 
as a function defined in th® £-plane satisfies the equation 


_d_ 

as 




(50) 


To determine the boundary conditions satisfied by cp 
observe that by virtue of the conformality of the mapping 
(28) with respect to the' metric (34) a line element normal 
to the profile P is taken into a line element normal to the 
circle | £ | = 1. Hence condition (20) implies that in the 

£-plane 


dcp 

dx 


<?(r 

cos T, r 

sin T) 

= 0 


dr 



r = i 


dcp a j + 

dcp Sn t 

dcp 

c is il + 

dcp dr» 

d | bx 

dn dx 

dy 

d % dy 

dn by 


( 51 ) 


Next 
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From (17) and (29) it follows- that 


t 




> 0 


lim 22-. 0 
£-> od &ri 




> 


J 


(52) 


From the symmetry of the flow (in fact, from condition (48)) 
it also follows that 


«p(| ,-n) = cp(| ,n) 


(53) 


Now let <Jj(£,Ti) denote the (one-valued) solution of eq.ua 
tion (50) satisfying the conditions 


_d_ 

dr 



r 


= 0 


1 


and 


lim 

08 


d$ 

dT 


1. 


lim 


d$ 

d-n 


0 


Then there exists a positive constant 0 


such that 


cp = 0 $ 


(54) 


(55) 


(56) 


This follows from the linearity and homogeneity of the differ- 
ential equations and the boundary conditions. 

The existence of the function $ is assured by known 
theorems and in principle $ can be computed numerically 
with any desired degree of accuracy. The actual computation 
of the function $ however, presents considerable diffi- 
culties. A brief discussion of the various possibilities 
will be found in appendix I. For the method described in 
this report it is irrelevant which numerical or mechanical 
device is actually used for solving equation (50), 
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Thus it is seen that the knowledge of the function f(u)) 
permits the computation of the potential cp as a function 
of £, except for a constant positive factor. 

Furthermore, to ever y function a = f(u>) satisfying 
conditions (31) and (32) there corresponds a function 
satisfying the condition 

$ = $(£,n) (5?) 

namely, the solution of (50) under the conditions (54) and 
(56) where T = T(|w*j 3 ) and w* = q*© -5 -® is the analytic 
function associated with f(uj) by means of equations (39) to 
(45). The function c£> can he computed if fw *** j nowhere 
exceeds unity and is determined by the shape function 0 (a) 
and by either of the two parameters q * or q* The 

actual determination of $ requires the numerical Integra- 
tion of a linear partial differential equation in a domain 
independent of the particular problem. 


6. The Functional Equation Satisfied by f(u>) 

Along the profile P the potential and the speed 

q are connected by the relation 


dcp 

q = a 0 — 
ds 

or 

ds = ^2. dcp 

q. 

If the function $ defined in the previous section is intro- 
duced, this can be written in the form 

fit c » 

= — - 2 — — -a_ $ ( COe UJ | S i n a,) du; (53) 

q Lf (uj ) J duu 

Setting 

= $(<ju) 


$(cos u), sin uj) 


(59) 
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integration yields 


I *- 


s = a 0 0 


$ * (u) 1 ) du) ' 

q[f(w‘)3 


But s = aS/2iT, q = $(q* 3 ), and cr = f(uj). 


Therefore 




f(ui) . so r 

2tt -J Q(OM) 


where 


$(w) = Q, I? 1 " [f(w)] S | 


By setting cu a 2 tt it is seen that 


SC 

2tt 


2tt 



. 5 2TT 


d$(o> ) 


so that finally 


f(ih) = 2 tt 


^.(JU dKoj* ) 
0 Q(o), 1 ) 

J.2TT d$(a) 1 ) 

0 <S(lO l ) 


(60) 


Since s.tarting from the function f(th) it is possible 
to compute q*(£) and therefore also q* [f(u>)3, Q,(u;), 

T(|,ri) and $>, equation (60) is a functional equation satis- 
fied hy the function f(<Ju) 
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7, Equivalence of the Functional Equation 
and the Poundary Value Problem 

In this section it will he shown that a solution of the 
functional equation actually yields a solution of the bound- 
ary value problem formulated in section 2. In the same time 
it will be shown how the knowledge of the function f(uu) 
permits the computation of the velocity distribution at 
points not situated on the profile P. 

In the following it is assumed that a function f (w ) 
satisfying (60) is known. This already implies that f is 
an increasing function (since its derivative is necessarily 
non*negative) and that the analytic function w * = q* e~^® 
associated with f(uu) (of. sec. 4) satisfies the condition 


| w*| < 1 (61) 

The function f (in ) determines functions q/a 0 , M, p/p 0 , 
T, 0 as functions of £ and Set 



where $ is the function defined in section 5, subscripts 
denote partial differentiation) and 0 is a positive constant. 
These line integrals are Independent of the path. (For the 
proof, see appendix II, sec. A), Hence eauation (62) defines 
a transformation 
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x = x(| ,Ti } , 


y = y (£ .*) 


(63) 


This transformation can also "be written in the form 


z » x + iy 


« Cl £o. e i6 - iO / . . a o . a 16 (ft- d|-$t dri) 

J * *^0/1” "M= £ 


(64) 


In order to compute the image of the circle |£| = 1 under 
this transformation the integration may he performed along 
this circle. But an the circle the normal derivative of $ 
Vanishes so that d £ ~ $£ dn - 0 and 


9 = 


( e[ftuj)] - TT for 0 < (JL) < TT 


1 


©[f(w)3 - 2 tt for tt < u> < tt 


^ = 5(u») 
q. 


i 

r- 

|5'| 

din 

for 

0 < in < it 

d$ = J 

1 

I IV . 1 




1 

L + 

1*1 

din 

for 

tt < in < 3tt 

so that the image of 

III 

= 1 

i s the 

curve 


z = C / e 
" o 


i0 [f (in' )] 1 d$ (w T ) I 
§(“■>') 


0 < in < 2tt 


or, hy (60) 


.in 


= °i / e i0 f i (y,! ) dcu> 0 < in ^ 


2tt 
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where 0, is a new positive constant. Thus, setting 
a = f((u) 


z = 



0 £ a < 2 


(65) 


Comparison of this with (12) shows that (62) takes the circle 
J£| = 1 into the profile P, and, since f(uo) is an in- 

creasing function, the mapping is one-to-one, The constant 
C x can he determined from the length of P, since this 
length is equal to 


S = C j 2tt 


(66) 


Next, the Jacohian of the transformation (62) equals 


d (x ,y ) 

d ( | ,n ) 



i 3 Jl 


1 


M a 




(67) 


It can he shown that the gradient of $ does not vanish for 
|£| >1. Hence the Jacohian is positive for |£| >1. Fi- 
nally, hy virtue of (55) and (46) as £ — * 


— 0 _ 0. 0 . — v C — 1 

dt qoo ’ 2)1 ’ dn * 1 „ - Mco 3 


where is the value of M at i = co . Hence 


( 68 ) 


lim z = cd (69) 

t— 


From the foregoing statements it follows that the trans- 
formation (62) yields a one-to-one mapping of the domain 
| £| > 1 into the domain E(P). Hence all functions defined 

in the £~plane q/a 0l 8, $, and so forth, can he considered 
as functions of x and y. The transformation inverse to 
( 63 ) 


i = S(x,y) , 


n = n (x , y ) 
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is conformal with, respect to the metric (34), The proof of 
this assertion will he found in appendix II, section B. 
Since $ satisfies equation (50) it follows from a lemma 
proved in a previous paper (reference 2; appendix, converse 
to Lemma l) that in the z-plane $ satisfies equat i on ” (15 ) 
considered as a linear differential equation (since p/p 0 

is a known function in the (-plane and therefore also in 
the z-plane). 

However, the relations 



a® 

+ 


by 

at 1 

dx 

at 

by 

bt 

d$ 

a® 

+ 


by 

Sri 

dx 

br\ 

by 

an 


together with (62) yield the equations 


$£ = 0 ^2. [[ cos $, 


4 l L 


Jl - M s 


*n 


$ + 
X 


sin 6$. - ° $^1 $ 

1 J 7 


s ■ ° v fr°° a e v t=m? ®s h + r*‘» + #=? ®i 


a 

q v 


L L 


*n 


/i-M= 


L 


yi - m : 


$ 


Solving for $ and $ yields 

A 

$ = ~ cos 0 , 

\J Q, 


Thus 


® y = 7T 7 1 - sin 9 
0 a 0 


u = q cos 8, v = q sin 8 

are the components of the gradient of the function 


ao'Pte.y) = 0 a 0 $ 
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The function cp satisfies equation (15) considered as a non- 
linear equation, that is, it satisfies equation ( 1 6 ) . Thus 
cp is the potential of a compressible flow, in the z-plane and 
u - iv = q e~*8 its complex velocity. 

It is now easy to verify that cp also satisfies the 
boundary conditions, In fact, along P the angle 9 coin- 
cides with the slope of the profile and at infinity 0 van- 
ishes. Finally, q takes the value q 0 = & 0 ^Cq*^) at 
infinity and its maximum is q maX = a Q Q, C q. max ^ * 

The foregoing discussion contains both the proof of the 
equivalence between the boundary value problem and the func- 
tional equation (60), and the description of a method permit- 
ting the computation of the flow in the whole domain B(P) 
once the function f (uu ) is known. 


8, Solution of the Functional Equation 

In the preceding seotions the determination of the pres- 
sure distribution along a given symmetrical profile P has 
been reduced to the determination of a function a - f(uu) 
satisfying the conditions 


f(0) a 0, f(2TT) = 2n 

f ( 2tt a) = 2 it - f (a ) 


(70) 


and the functional equation (60). The form of this equation 
suggests the application of the iteration method. 


Let f(cn) be any function satisfying conditions (70). 

By use of this function the functions w*, §(u>) , T and §(uu) 
can be computed (cf. secs. 4 and 5) provided the function w* 
satisfies the condition I w*| < 1, (This will always be the 

case if is prescribed.) Therefore the right-hand side 

of (60) can be computed and will represent an increasing 
function g(uu) which will again satisfy conditions (70), 

This function will be denoted by 

g(w) = F [f (u>‘ ), u>] (71) 

With this notation equation (60) reads 

f (tu) = F [f (uj 1 ) , u>] 


(72) 
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Thus the problem consists of finding a fixed point of the 
functional transformation IP. 

Now let f 0 (uu) (approximation of order zero) be a 
function satisfying conditions (70) and Bet 


f i (u>) = F [f 0 (uj * ) , uj] 
f 3 (u>) = F [f x (o)» ) , uj] 




> (73) 


f n (uj) = F lf n _ 1 (0J> ) , uj] 




If the sequence f Q (uj), f x (cu ) , , , ,, f n (ui), » . , converges 

to a function f(uj) and if lim F( n ) = F(f), then f(uj) is 
a solution of equation (60). _ 

It seems to be rather difficult to prove rigorously the 
convergence of this method. However, on the basis of computa- 
tions carried out for the special case Y * -1 a fairly 
rapid convergence of the iterations seems probable, and it 
might be expected that good results will be obtained after a 
few steps. The real difficulty in computing the successive 
approximations lies in the numerical integration of the dif- 
ferential equation satisfied by <£. (Of. appendix I.) Thus 
the practical applicability of this method hinges upon the 
development of efficient methods for solving this problem, 
which, of course, is of prime importance for practically all 
branches of engineering mathematics. . 

It should be noted that the iteration process may break 
down if in the original boundary value problem the value of 
<3.0, is prescribed. This cannot happen if the value of q ma3t; 

is prescribed. 

The rapidity of the convergence will depend upon the 
appropriate choice of the function f 0 (u>) (the approxima- 
tion of order zero). Several possibilities of choosing f 0 
are listed in order of preference. 
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(1) Choose for £ 0 the solution of equation (60) pre- 
viously obtained, for a profile P* close to the desired pro- 
file and for a value of q’ (or q* m«x) close to the desired 
value of q tt (or q aaX ) . 

(2) Choose for f 0 the solution of equation (60) for 
the special case Y = -1. (Cf. next sec. ) 

(3) Choose for fo the function resulting from the con- 
formal mapping of the profile into a circle. 

(4) Por thin profiles £ 0 (u>) may be chosen as 


f 0 (u>) = j (1 - cos U)) , 
f Q (2TT - uj) = 2tt - f (tjj) 

J 

Note that (3) and (4) are special cases of ( 1 ). (in the 

case (3) P ! *e P and q ^ - 0, in the case (4) P' is a 
straight segment and q ^ = 0.) 


0 < O) < TT 




(74) 


9, Comparison with the Case Y « -1 

If it is assumed that the fluid obeys the so-called lin- 
earized pressure-density relation of Chaplygin-Karman-Ts i en 
(references 4, 6, and 7), Y = -1 and equations (l), (3), 

and (4) yield 


T = 1 





) 


It is seen that in this case the function $ becomes inde- 
pendent of the particular function f(ui). In fact, 


0 = 1 + 


e 

I s + -n 3 


so that 
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$ (uj) = 2 COS 


UJ 


Equation (60) takes the form of an integral equation treated 
in a previous paper (reference 1, equation (5l)). Since 
does not have to he determined, only functions of one vari- 
able enter in the problem and the computational labor is 
greatly diminished. The convergence of the iteration method 
described in the preceding section when applied to this inte- 
gral equation turned out to be very satisfactory. 


10, Comparison with Other Methods 

All other analytical methods proposed for solving the 
boundary value problem of section 2 involve an expansion of 
the potential function cp(x,y) in terms of a parameter char- 
acterizing either the deviation of the flow from the uniform 
flow (such as the thicxness parameter of a profile) or the 
deviation of the flow from an incompressible flow (such ae 
the stream Mach number). These methods also involve succes- 
sive integrations of linear partial differential equations. 
Methods based on such expansions necessarily suffer from the 
disadvantage that the convergence becomes worse as the devi- 
ation from the "undisturbed state" increases, that is, pre- 
cisely as the influence of compressibility becomes more pro- 
nounced. In the present method the incompressible (i.e., 
infinitely slow) flow or the uniform flow (flow past a 
straight segment) is in no way distinguished, and the rate of 
convergence should not depend too muoh upon the value of the 
pa r s ^ ^max* 

Furthermore, in all other methods the speed at infinity 
and therefore the stream Mach number are prescribed, and 
there is no way of telling whether this value of the stream 
Mach number does not exceed the critical Mach number. If 
this is the case, the flow becomes partly supersonic, and 
the convergence usually breaks down. In the present method 
it is possible to prescribe the maximum local Mach number 
M max and the value of M maz may be chosen arbitrarily close 
to unity. 

The preceding remark applies also to a purely numerical 
treatment of the nonlinear equation (15) (see a recent paper 
by Salmons , reference 8, and the literature quoted therein;. 
This treatment consists in replacing the differential equa- 
tion by a nonlinear difference equation which is solved by 
the relaxation method. 
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11. On a Velocity Correction Formula by Garriok and. Kaplan 

The functional equation to which the pressure distribu- 
tion problem has been reduced throws new light on a velocity 
correction formula recently proposed by G-arrick and Kaplan. 
(See reference 9.) These authors derived their "geometric 
mean formula" using an analogy between compressible and in- 
compressible flows and some ideas from the theory of sigma- 
monogenic functions. In order to write this correction for- 
mula in the notations of the present paper, let XCM^) denote 

the value of q* corresponding to the stream Mach number 
and set 

<^(q*) = Q(<f a ) (76) 

Then the Garrick-Kaplan formula reads 




[* <M J ] 


(76) 


Here (q/q w ) denotes the ratio of the speed at a point of 
a profile to the speed at infinity for a compressible flow 
of stream Mach number and ^ ( l/ c loo^i is the same ratio 

'for. aa ..incompressible flow (i.e,, for a flow with = 0). 

For the case of a fluid obeying the linearized pressure- 
density relation with Y — -1 formula (76) takes the form 

/'~ s \ » f. JL' \ } , ~ .,. h ( 77 ) 

W 0 Vi 1 - ^(q/qji 

where 

* - ^(M w ) = — (78) 

i +yi - m a 

v a> 

This is the Karman-Tsien velocity correction formula. (See 
reference 6. ) 

It is immediately seen that formula (76) would be rigor- 
ously correct if the function a = f 0 (u>) resulting from the 
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conformal mapping of the profile into a circle were the solu- 
tion of the functional equation (60) for the desired value 
of M^, The situation may he summed up as follows! The dif- 
ference between the relative velocity distribution (the dis- 
tribution of the values q/q^) in an incompressible flow 
and in a compressible flow of stream Mach number is due 

(a) to the difference between the values of the corresponding 
mapping functions a = f (u> ) , and (b) to the fact that the way 
in which the velocity distribution is determined by the map- 
ping function f(ui) depends upon the stream Mach number. 

The velocity correction formula (76) takes into account the 
factor (b) but fails to include the effect of factor (a). 

On the basis of computations carried out for the case 
y » -1 (see reference 1, figs. 1 and 2) it may be concluded 
that formula (76) will yield values of q/q^ which are too 

low. Therefore the values of the critical Mach number 
(stream Mach number for which = l) computed by formula 

(76) may be expected to be higher than the actual values. 


CONCLUDING REMARKS 


1, The method presented in this report reduces the 
problem of finding the velocity and pressure distribution 
along a preassigned symmetrical profile to the solution of a 
certain mapping problem which in turn is reduced to the solu- 
tion of a functional equation by the method of successive ap- 
proximations. Though the mathematical difficulties peculiar 
to the nonlinear differential equations of fluid dynamics are 
eliminated by this method, the practical application of this 
procedure requires the numerical solution of linear partial 
differential equations with variable coefficients. This 
last problem presents no theoretical difficulties, but really 
efficient methods for its solution are still in the process 
of development, Upon this development, which is intimately 
connected with the design of automatic computing machines, 
hinges the practical applicability of the present method. 

2, The extension of the present method to the case of 
circulatory flows and to the case of flows with locally 
supersonic regions will require a more profound investiga- 
tion of the "mappings conformal with respect to a compress- 
ible flow" introduced by the author in a previous report. 

(See reference 3,) 
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3. The author take* the opportunity to draw attention 
to the similarity between the method given in reference 3 
and the one given (from a different point of view) in a 
paper by S. A. Chrietianovitoh (reference 17). In the 
author's opinion, however, the paper "by the noted Russian 
fluid dynamicist contains an error. In the terminology 
introduced in reference S y a statement made by Chrietianovitoh 
oould be formulated as follows! "To every subsonic compress- 
ible flow past a profile P there exists a conjugate (modulo 
l) incompressible flow pas-t another profile P 1 ." The author 
believes that this statement i* misleading and that, for this 
reason, Christianovitch** method muet fail in the case of 
circulatory flows. 


Brown University, 

Providence, R. I., S*ptenifc*r 1&&B . 


iPPMTS iit t 

CONCERNING- THE COMPUTATION OP THE FUNCTION $ 


A, No matter which method is used for integrating the’ 
linear partial differential equation (50) satisfied by $ 
it is convenient first to introduce the auxiliary complex 
variable Z = X + iY by the relation 

Z = t + ^ (Al) 

This transformation takes thB domain ] £l > 1 into the do- 
main exterior to the straight segment (-1, 4-1 ) . The func- 
tions $ and T may be considered as functions of X and 
Y, and in view of the conformality of the transformation 
(Al) $ satisfies the equation 


- ( 

1 


(k 

dx \ 

.T aX> 

&Y 

\T 5Y/ 


Furthermore, by virtue of (57) 


(A3) 


$ (X, ~Y) * $ (X , Y ) 


(A3) 
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This implies that 


Y = 0 


(A4) 


for X < -1 and. X > 1. But the same relation holds for 
-1 < X < 1 by virtue of (54). At infinity equation (55) 
implies that 


1 M_^o as Z — * » 

ax 6Y 


( A5 ) 


The determination of 4> requires the integration of equa- 
tion ( A2 ) in the half-plane Y > 0 under the boundary con- 
ditions (A4) and (A5 ) , 

A numerical or mechanical integration can he performed 
in a finite domain only. Therefore the upper half-plane is 
replaced by a sufficiently large rectangle with the vertices 
(-A, 0), (-A, B), (A, B), (A, 0) and the condition (A5) is 

replaced by the boundary conditions 

$ (±A, Y) » ±A, $ (X, B) = X (A6) 

(Another way of passing to a finite domain would consist of 
a conformal transformation, say of |£| > 1 into | £j < 1.) 

B, The boundary value stated formerly can be solved by 
means of electrical models. One method is that of the elec- 
trolytic bath due to Taylor and Sharman (references 10 and 11), 
the other that of network analogies due to Kron (references 12 
and 13). 

0, In order to solve equation (A3) numerically it is 
first replaced by a difference equation. (The possibility 
of reducing (A2) to an integral equation will not be dis- 
cussed here.) Let X 0 , Y 0 be a point of the domain con- 
sidered and set X a = X 0 + 8 , Yj = Y 0 , X 3 = X 0 , Y s = Y 0 + fi, 

^3 ~ Y o ** 8 t Y 3 = Y o* ^4 = i Y 4 = Y o “ 8 , 8 being a small 
positive number. Denote the value of a function Q at 
( x i» 7 ±) "by and replace the expressions 
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by 

ci 1 - n 3 a, - a 4 + n 3 - 2 n 0 n 3 + n 4 - 2Q 0 

f , ~ > s 

28 28 8 2 * 6 a 

respectively, Equation (A2) can be rewritten in the form 

4T 0 ($ 1 + $ 3 + $ 3 + 4$ 0 ) = (T x - T 3 )(cp 1 ~ tp 3 )+ (T a - T 4 ) (cp a - cp 4 ) 


or 


cp 0 = e lCPl + 3 2 cp a + 3 3 q> : 


( A7 ) 


where 


Pi = 


3, = 


4T 0 - T 1 + T 3 
16 T„ 


4To + Ti - T 3 


16 T, 


3 2 = 


3. = 


4T 0 - T g + 


16 T, 


4T + T - T a 
o a 4 


\ (A8) 


16 T 


o J 


If the point 0 lies on the line Y = 0 , then the point 
4 lies outside of the domain considered and according to 
(A3) the preceding formulas must be changed to 


\ (A9) 

P-i - 0 

J 

Note that the following relation always holds 

p i + P a + P 3 + P 4 = 1 

and that 

£ 0 (A10) 


Pi = 


4T q - + T 3 

16 T n 


P, = 


4T 0 + T x - T ; 
16 T 


if 8 is sufficiently small 



NACA TN No. 1012 


36 


Since equation (A2) is invariant with respect to rota- 
tion, the same difference equation would he obtained if the 
square formed by the points 1, 2, 3, and 4 would not have 
its sides parallel to the coordinate axes. Therefore, it 
is easy to compute $ at lattice points which do not form a 
uniform net. In fact, it is advisable to use more points 
near and on the segment (-1, l) for here the function T 
will change rapidly and the values of $ on this segment 
are the ones actually needed. 

A simple computation shows that if the function 


* (X, Y) . Y c nn x“ Y m 


satisfies the differential equation (a 2), then the polynomial 

®oo + °io^ + c oi ^ + °8 0 ^ 3 + + °oa^ 3 

satisfies the difference equation (A7). It also is known 
that solutions of the difference equation converge toward 
the solution of the differential equation as the parameter 
6 approaches zero. 

D, The boundary value problem for the difference equa- 
tion derived above can be solved either by a rigid iteration 
scheme, say by Liebmann's method (references 14 and 15) or 
by a relaxation method, Concerning the latter, reference Is 
made to a recent paper by Emmons (reference 8) and to the 
literature quoted therein. The iteration methods require 
more computational labor than the consequent application of 
the relaxation method and the use of modern computing devices 
becomes imperative. On the other hand, the instructions for 
computing become simpler and the whole problem can be "coded" 
once and for all on an automatic computing machine. An at- 
tempt was made to solve equation (A2) by Liebmann’s method 
using the I.B.M, Automatic Sequence Controlled Calculator, 

This attempt could not be considered successful because of 
the slowness of convergence. It is possible that convergence 
could be improved by using either a difference equation in- 
volving higher differences or a different Iteration procedure. 
However, the design of automatic computing machines is now in 
the process of rapid development, and the speed at which 
these machines perform will doubtless be increased at a 
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tremendous rate, This trend will, on the one hand, reduce 
the importance of the rate of conversance and, on the other, 
put a premium on the simplicity of the computational proced- 
ure. 


E. In this section, the convergence of Liebmann's itera- 
tion method when applied to the boundary value problem for 
the difference equations (A7 ) to (A3^ will be proved, The 
proof is patterned after the one given by Loewner for the 
case of the first boundary value problem for Laplace’s dif- 
ference equation. (See reference 15,) 


Let P lt P a , . . . , Pjj be the points at which the val- 
ues of $ (denoted by $ a , . . ., are to be computed 

It will be assumed that among the four neighboring points of 
the point P* there Is at least one boundary point (i.e., a 
point with either x = ^A or Y « B), and that one of the 
four neighboring points of P n is P^^n * 2 . 3 . • • •. N ). 

Liehmann’s procedure for finding the values is as fol- 

lows. Assume some trial values Oj 0 , 

Go over the points P^ in the indicated order correcting the 
value <5^°^ using equations (A?) and (A 8 ) (or equations 

(A7 ) and (A9) if P^ lies on the X-axis). In applying these 
formulas use corrected values of $ whenever possible. In 


this way new values $ 3 ^, .... are obtained 

(a ) (s ) 

Repeat this process to obtain the values $4 1 , , ... 

It iB asserted that as n — *■ 00 11 converges toward the 

desired value 


In what follows it is assumed that the net is so fine 
that condition (A10) is satisfied. Plainly 


and , 


*i 


U) 


in general. 


X 4 J <0) + 


( a+ 1 ) 
®i 


I 


(n) 


(All) 


p ij $ i 


+ ^io 


(A12) 
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where the constants 3 ^ . are linear combinations of products 
of the 3 t defined by (A 8 ) and (A9) and the constants 3 i0 
depend upon the boundary conditions. Thus 

j > 0 ( A13 ) 

Now assume that the boundary conditions read 

$ * 0 for X = ^A and for Y = B 

and that = 1 , i ■ 1 , 2 , . . . N. It is easily seen 

that in this case all 3 * „ vanish and that < $i • 

* (») ^ „ (o) 10 

®a < , . . . . Since, in the case considered 

^ <l) ■ x ^ 

it follows that there exists a constant b such that 

3 1<} < b < 1, i = 1, 2, . . . , N (A14) 

J 

Since, in the general case 

$ i (n) = ^ 8 i<3 + ki (A16) 

then from (A12) and (A15) it follows that 

$ i (n+1 > - Pij (V n) - (A 16 ) 

Set 

. - (n) . (n~i) 

d n = max - $ 3 . CA17; 

Jroa (A16) and (A13) it follows that 
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r I 


so that by (A14) 


d n+x * d n maX 


d n+i < * d n 


U 


and hence 


n-i 

<i n < * d i 


as n > co 


(A18) 


Thus the limits 


<$. = lim 
1 n — * 05 i 


(n) 


exist, since 


(A19) 


(n+m) 

(n) 


(n+m-p) 

(n+m-p-i ) 

<6 

v i 

" ®i 

Z, 

p=0 

*i 

“ ®i 


m-x 

< 1 * 

p=o 


n+m-p-x 


dj, * d a b 


n i „ -b m d x b 


1 - h 1-1) 


Plainly 




t ■ Z 


J + k i 


(A20) 


which shows that the 0^ satisfy both the boundary condi- 
tions and the difference equation. 


Also, from (A12) and (A20) it follows that 

4 i ( ° +1> - *i - I “ij (®J <n) " *4 


) 


so that by (A13 ) and (A14) 
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max 

i 



(n+i ) 




< max 
i 



( A21 } 


This remark establishes the uniqueness of the solution. 

Finally, if all differences ar« posi- 

tive (negative), then by (A16) and (A13) so are all the dif- 
ferences $ ^ n+1 ^ This implies that all values 

lie below (above) the desired values 

To summarise, the following statements are seen to be 

true . 


(i) Liebmann's method converges, 

(ii) Both the maximum corrections 

(n) 


and the maximum differences 
ically . 


$ 


(n+i ) 


- 


(n) 


(iii) If all corrections $ 


®i - ®i 

(n+i ) 


i " v i 

decrease monoton- 


— $. 


(n) 


are positive 


/ V 1 A 

(negative) all values n+1 are smaller (larger) than the 

desired values. 


If no assumptions concerning the order of the points 
are made, a slight refinement of the argument shows that 
statements (i) and (iii) remain true and (ii) holds if the 
words "decrease monotoni cally " are replaced by "do not in- 
oreas e . " 


APPENDIX II 


CONCERNING THE MAPPING FROM THE Z-PLANE INTO THE £- PLANE 


A. This appendix contains the proof of two assertions 
concerning the transformation (62) made in section 7, First, 
the independence of the line integrals on the right-hand side 
of (62) on the path of Integration will be established. Sine 
$ satisfies equation (60) there exists a function 'K£,n) 
such that 
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= T 

^ = -I V|/g 


( AS2 ) 


(The function is proportional to the stream-function 

of the compressible flow in the £ -plane.) Since log w* is 
an analytic function of £, $ and may "be considered as 

functions of log q* and 8. Since equations (A22) are in- 
variant with respect to conformal transformations, it fol- 
lows that 


5 8 5 1 og q* 

M = _ r 

5 log q* 5 8 


(A23) 


By introducing instead of log q* the new variable q con- 
nected with log q* by equation (3) and using equation (4), 
the following system is obtained 

5$ p 0 q 5\|/_ 

58 p 5q 

( A24) 

5£ a p 0 (i - M a ) 

5q pq 56 


(These are the well-known Chaplygin equat 
Kow (62) may be written in the form 


I 3 0 


C 0 8 


Po 

p 


sin 8ty|^ d| + (^eoi 


y a o J"* ~~ ~ {^C sin + 008 0^1^ 4 | + ^si] 


ions . ) 


6$ r ! -~- 8in 0^r| 
0$ + ^2- cos 

'n p *n 


>”} 

>"} 


or 
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lx 


dy 



0 i { ■ sin 

8 4 $ - cos 


e d \|^ 

6d 




( A25 ) 


and equations (A24) are precisely the conditions that these 
integrals he independent of the path (cf. reference 16, p. 
8 ). 


B. It will now he shown that the mapping inverse to 
(62) is conformal with respect to the metric (34), (35). By 
use of the relations 


dz 

dh 

T 

dx 

jiL 

H x 

dy ' 

dri 

“ J dy 

dy 

- j in, 

dy 

. 

ae 

dx' 

' dn 

dx 


j - d (* ,y) 
d(£,n) 


together with (62) and (67), the following equations are ob- 
tained 


« ($|+ #®)C/l - M 8 sin 6®^ + cos 0®|) 

( + (-v / l"" - M 8 cos 6$r| + sin6$| ) 


dz 

ii 

dy 

dr> 


— = (4>| + 4^) {-Jl - M S sin 0$£ + cos ©4^) 


= (4>t + ®*) {>Jl - K 3 cos 04». + sin 04> ) 


dy 


'ft 


7) 
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It is seen that these partial derivatives satisfy the condi- 
tions 


dn 



where the g^ fc are given hy (35). These are the well-known 
Beltrami equations expressing the conformality of the mapping 

t 3 Kx,y ) , n * n(x,y ) 

with respect to the metric (34). 
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. 3 -1 \ f, . r . . 4 


TABLE l. r ’ q*2 and T as functions of (l-li 8 



I 


.00 

.01 

.02 

0.03 

0.04 


0.05 

0.06 

0.07 

0.08 

0.09 

0.10 

0.11 

0.12 

0.15 

0.14 

0.15 
.16 
.17 
.18 
0.19 

0.20 

0.21 

0.22 

0.25 

0.24 

0.25 

0.26 

0.27 

0.28 

0,29 

0.50 

0.51 

0.52 

0.55 


1.00000. 
1.00000 
.U 00000 
.-0.99999 
0.99996 

Q. 99995 
: 0.99988 
0. 99981 
0.99971 
• 0.99959 

0.99944 

0.99926 

0.99905 

0.99877 

0.99846 

1 6.99810 
0.99769 
0.99722 
0.99670 
0.99611 

0.99545 

0.99472 

0.99591 

0.99502 

0.99204 

0.99098 

0.98982 

0.98856 

0.98720 

0.98575 

0.98415 

0.98244 

0.98061 

0.97865 


0.00000 
0.0lfc77 
0 . Q5155 
0.04728 
0.06502 

0.07875 

0.09446 

0.11014 

0.12580 

0.14142 

0.15701 

0.17256 

0.18807 

0.20355 

0.21894 

0.23429 

0.24959 

0.26482 

0.27999 

0.29509 

0.31012 

0.32506 

0.33993 

0.35471 

0.36940 

0.38400 

0.39851 

0.41891 

0.42722 

0.44141 

0.45550 

0.46947 

0.48332 

0.49706 



0.34 

0.35- 

0.36 

0.37 

0.38 

0.30 

0.40 

0.41 

0.42 

0.43 

0.44 

0.46 

0.46 

0,47 

0.48 

0.49 

0,50 

0.51 

0.52 

0.55 

0,54 

0.55 

0.56 

0.57 

0.58 

0.59 

0.60 

0.61 

0.62 

0.63 

0.64 

0.65 

0.66 

0.67 


0.97656 

0.97432 

0.97194 

0.96939 

0.96670 

0.96583 

0.96079 

0.95758 

0.95417 

0.95057 

0.94676 

0.94275 

0.93852 

0.93406 

0.92936 

0. 92443 
•0.91923 
O; 91378 
0.90806 
0.90205 

0.89575 

0.88915 

0.88223 

0.87499 

0.86741 

0.85948 

0.85119 

0.84252 

0.83346 

0.82401 

0.81413 

0.80382 

0.79306 

0.78184 
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TABLE 1. (Continued) 



1-M 2 

(9*) 8 

T 

0..68 

0.77014 

0.87765 

0.69 

0.75794 

0.88506 

0.70 

0.74523 

0.89226 

0.71 

0.73198 

0.89923 

0.72 

0.71817 

0.90598 

0 .73 

0.70380 

0.91251 

0.74 

0.68882 

0.91881 

0.75 

0.67324 

0.92489 

0.76 

0.65701 

0.93073 

0.77 

0.64012 

0.93635 

0.78 

0.62255 

0.94173 

0.79 

0.60427 

0.94688 

0.80 

0.58525 

0.95180 

0.81 

0.56546 

0.95648 

0.82 

0.54490 

0.96093 

0.83 

0.52351 

0.96514 

0.84 

0.50127 

0.96911 

0.85 

0.47814 

0.97285 

0.86 

0.45411 

0.97634 

0.87 

0.42913 

0.97960 

0.88 

0.40316 

0.98262 

0.89 

0.37618 • 

0.98540 

0.90 

0.- 34813 

0,98793 

0.91 

0.31900 

0.99023 

0.92 

0.28872 

0.99228 

0.93 

0.25726 

0.99409 

0.94 

0.22458 

0.99566 

0.95 

0.19062 

0.99699 

0.96 

0.15535 

0.99807 

0.97 

0.11870 

0.99892 

0.98 

0.08063 

0.99952 

0.99 

0.04108 ' 

0.99988 

1.00 

0.00000 

1.00000 
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TABLE 2. 1-M S and T as functions of q*2 


( q *) 2 

0.000 

0.005 

0.010 

0.015 

0.020 

0.025 

0.030 

0.035 

0.040 

0.045 

0.050 
0.055 
060 
065 
070 

075 
080 
085 
0.090 
0.095 

0.100 

0.105 

0.110 

0.115 

0.120 

0.125 

0.130 

0.135 

0.140 

0.145 

0.150 

0.155 

0.160 

0.165 


\ 1-M 2 

T 

(q*) 2 

■esi 

1.00000 

’ 1.00000 

0.170 

0.95589 

0.99880 

1.00000 

0.175 

0.95448 

0.99760 

0.99999 

0.180 

0.95305 

0.99639 

0.99998 

0.-185 

0.95162 

0.99518 

0.99997 

0.190 

0.95018 

0.99396 

0.99996 

0.195 

0.94873 

0.99274 

0.99994 

0.200 

0.94728 

0.99151 

0.99991 

0,205 

0.94581 

0.99027 

0.99989 

0.210 

0.94434 

0.98903 

0.99986 

0.215 

0.94286 

0.98778 

0.99982 

0.220 

0.94137 

0.98652 

0.99978 

0.225 

0.93987 

0.98526 

0.99974 

0.230 

0.93837 

0.98400 

0.99969 

0.235 

0.93685 

0.98273 

0.99964 

0.240 

0.93533 

0.98145 

0.99959 

0.245 

0.93380 

0.98016 

0.99953 

O'. 250 

0.93226 

0.97887 

0,99946 

0.255 

0.93070 

0.97757 

0,99939 

0.260 

0.92914 

0. 97627 

0.99932 

0.265 

0.92758 

0.97496 

0.99924 

0.270 

0.92600 

0.97364 

0.99916 

0.875 

0.92441 

0 . 97232 

0.99908 

0.280 

0.92281 

0.97099 

0.99899 

0.285 

0.92120 

0.96965 

0.99889 

0.290 

0.91958 

0.96831 

0.99879 

0.295 

0.91796 

0.96696 

0.99868 

0.300 

0.91632 

0.96560 

0.99857 

0.305 

0. 91467 

0.96424 

0.99846 

0.310 

0.91301 

0.96286 

0.99834 

0.315 

0.91134 

0.96148 

0.99821 

0.320 

0.90966 

0.96010 

0.99808 

0.325 

•0.90797 

0.95870 

0.99794 

0.330 

0.90627 

0.95730 

0.99780 

0.335 

0.340 

0.90455 

0.90283 


T 


0.99766 

0.99750 

0.99734 

0.99718 

0.99701 

0.99683 

0.99665 

0.99646 

0.99626 

0.99606 

0.99586 
0.99564 
0. 99542 
0.99519 
0. 99496 

0.99471 

0,99446 

0.99421 

0.99394 

0.99367 

0.99339 
0.99311 
0.99281 
0.99251 
0. 99220 

0.99188 

0.99155 

0.99121 

0.99087 

0.99051 

0.99015 

0.98978 

0.98940 

0.98900 

0.98860 
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TABLE 2 . ( C ont inued) . 


(q*) 2 

HUB 

T 


(q*) 2 

■59 

T 

0.345 

0.90109 

0.98819- 


0.500 

0.84056 

0.96933 

0.350 

0.89935 

0.98777 


0.505 

0.83835 

0.96847 

0.355 

0.89759 

0.98734 


0.510 

0.83612 

0.96760 

0.360 

0.89582 

0.98690 


0.515 

0.83387 

0.96671 

0.365 

0.89403 

0.98645 


0.520 

0.83160 

0.96579 

0.370 

0.89224 

0.98598 


0.525 

0.82931 

0.96486 

0.375 

0.89043 

0.98551 


0.530 

0.82701 

0.96390 

0.380 

0.88861 

0.98502 


0.535 

0.82468 

0.96293 

0.385 

0.88677 

0.98453 


0.540 

0.82232 

0.96193 

0.390 

0.88492 

0.98402 


0,545 

0.81995 

0.96091 

0.395 

0.88306 

0.98350 


0.550 

0.81756 

0.95986 

0.400 

0.88119 

0.98296 


0.555 

0.81514 

0.95880 

0.405 

0.87930 

0.98242 


0.560 

0.81270 

0.95770 

0.410 

0.87740 

0.98186 


0.565 

0.81023 

0.95659 

0.415 

0.87549 

0.98129 


0.570 

0.80774 

0.95545 

0.420 

0.87356 

0.98070 


0.575 

0.80523 

0.95428 

0.425 

0.87161 

0.98011 


0.580 

0.80269 

0.95308 

0.430 

0.86966 

0.97949 


0.585 

0.80013 

0.95186 

0.435 

0.86768 

0.97887 


0.590 

0.79754 

0.95061 

0.440 

0.86569 

0.97823 


0.595 

0.79492 

0.94933 

0.445 

0.86369 

0.97757 


0.600 

0.79228 

0.94802 

0.450 

0.86167 

0.97691 


0.605 

0.78961 

0.94668 

0.455 

0.85964 

0.97622 


0.610 

0.78691 

0.94531 

0.460 

0.85758 

0.97552 


0.615 

0.78418 

0.94391 

0.465 

0.85552 

0.97481 


0.620 

0.78142 

0.94247 

0.470 

0.85343 

0.97408 


0.625 

0.77863 

0.94101 

0.475 

0.85133 

0.97333 


0.630 

O'. 77581 

0.93950 

0.480 

0.84921 

0.97256 


0.635 

0.77296 

0.93796 

0.485 

0.84708 

0.97178 


0.640 

0.77007 

0.93639 

0.490 

0.84492 

0.97098 


0.645 

0.76715 

0.93477 

0.495 

0.84275 

0.97016 


0.650 

0.76420 

0.93312 
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TABLE 2 . {Continued) 


(q*) 2 

|l-M S 

T 

0.655 

0.76121 

• 0*93142 

0.660 

0.75819 

0*92969 

0.665 

0.75513 

0.92791 

0.670 

0.75203 

0.92609 

0.675 

0.74889 

0.92422 

0.680 

0.74571 

0.92231 

0.685 

0.74249 

0.92035 

0.690 

0.73923 

0.91833 

0.695 

0.73592 

0.91627 

0.700 

' 0.73257 

0.91415 

0.705 

0.72918 

0.91198 

0.710 

0.72573 

0.90975 

0.715 

0.72224 

0.90747 

0.720 

0.71870 

0.90512 

0.725 

0.71511 

0.9027.1 

0.730 

0.71146 

0.90023 

0.735 

0.70775 

0.89768 

0.740 

0.70399 

0.89507 

0.745 

0.70017 

0.89238 

0.750 

0.69629 

0.88962 

0.755 

0.69235 

0.88677 

0.760 

0.68834 

0.88385 

0.765 

0.68427 

0.88084 

0.770 

0.68012 

0.87774 

0.775 

0.67590 

0.87454 

0.780 

0*67160 

0.87125 

0.785 

6.66723 

0.86786 

0.790 

0.66277 

0.86436 

0.795 

0.65823 

0.86075 

0.800 

0.65360 

0.85703 

0.805 

0.64888 

0.85318 

0.810 

0.64406 

0.84921 

0.815 

0.63914 

0.84510 

0.820 

0.63411 

0.84085 

0.825 

0.62897 

0.83645 

0.830 

0.62371 

0.83190 

0.835 

0.61833 

0.82718 

0.840 

0.61283 

0.82228 

0.845 

0.60718 

0.81720 

0.850 

0.60140 

0.81193 


(q*) 2 

^1-M S 

T 

0.855 

0.59546 

0.80644 

0.860 

0.58936 

0.80073 

0.865 

0.58309 

0.79479 

0.870 

0.57663 

0 . 788/59 

0.875 

0.56998 

0.78212 

0.880 

0.56313 

0.77535 

0.885 

0.55605 

0.76827 

0.890 

0.54874 

0.76086 

0.895 

0.54116 

0.75307 

0.900 

0 . 53331 

0.74488 

0.905 

0.52515 

0.73626 

0.910 

0.51666 

0.72715 

0.915 

0.50781 

0.71752 

0.920 

0.49856 

0.70730 

0.925 

0.48886 

0.69643 

0.930 

0.47867 

0.68483 

0.935 

0.46793 

0.67241 

0.940 

0.45655 

0.65905 

0.945 

0.44446 j 

0.64460 

0.950 

0.43152 : 

0.62890 

0.955 

0.41761 

0.61171 

0.960 

0.40252 

0.59274 

0.965 

0.38600 

0. 5^159 

0.970 

0.36769 

0.54767 

0.975 

0.34704 

0.52017 

0.980 

0.32320 

0.48774 

0.985 

0.29471 

0.44806 

0.990 

0.25852 

0.39636 

0.995 

0.20623 

0.31944 

1.000 

0.00000 

0.00000 
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TABLE 5. M p/p 0 , and q/£ 0 as functions of q* s 


(q*>® 

M 

P ' s * 

p/p c 

q/a Q 

0.00 

0.0000 

1.0000 

1.0600 

0.0000 

0.02 

O'. 0982 

0.9952 

0.9933 

0.0981 

0.04 

0.1392 

0.9904 

0.9865 

0.1389 

0.06 

0.1711 

0.9855 

0.9797 

0.1706 

0.08 

0.1982 

0.9806 

0.9729 

0.1974 

0.10 

0.2224 

0.9757 

0.9660 

0.2213 

0.12 

0.2445 

0.9707 

0.9591 

0.2430 

0.14 

0.2650 

0.9657 

0.9522 

0.2631 

0.16 

0.2844 

0.9607 

0.9452 

0.2821 

0.18 

0.3028 

0.9556 

0.9382 

0.3000 

0.20 

0.3204 

0.9505 

0.9311 

0.3171 

.0.22 

0.3374 

0.9453 

0.9240 

0.3336 

0.24 

0.3538 

0.9401 

0.9168 

0.3494 

0.26 

0.3678 

0.9348 

0.9096 

0.3648 

0.28 

0.3852 

0.9295 

0.9024 

0.3795 

0.30 

0.4005 

0.9241 

0.8951 

0.3942 

0.32 

0.4154 

0.9187 

0.8878 

0.4083 

0.34 

0.4300 

0.9132 

0.8803 

0.4222 

0.36 

0.4444 

0.9077 

0.8728 

0.4358 

0.38 

0.4587 

0.9021 

0.8653 

0.4492 

0.40 

0.4728 

0.8965 

0.8576 

0.4625 

0.42 

0.4867 

0.8908 

0.8500 

0.4754 

0.44 

0.5005 

0.8850 

0.8422 

0.4883 

0.46 

0.5143 

0.8791 

0.8344 

0.5011 

0.48 

0.5281 

0.8732 

0.8265 

0.5138 

0.50 

0:541? 

• 0.8672 

0.8185 

0.5263 

0.52 

0.5554 

0.8611 

0.8104 

0.5388 

0.54* 

0.5690 

0.8549 

0.8022 

0.5512 

0.56 

0.5827 

0.8486 

0.7940 

0.5637 

0.58 

0.5964 

0.8422 

0.7856 

0.5760 

0.60 

0.6102 

0 . 8357 

0.7771 

0.5884 

0,62 

0.6240 

0.8291 

0.7685 

0.6008 

0.64 

0.6380 

. 0.8224 

0.7598 

0.6132 

0.66 

.0.6520 

0.8155 

0.7509 

0.6256 

0.68 

0.6663 

0.8085 

0.7418 

0.6382 
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TABLE 3 . (Continued) 


(< i *) 2 

M 

P f Po 

p / p 0 

« l / a 0 

0.70 

0.680 r / 

0.8014 

0.7326 

0.6508 

0.72 

0.6953 

0.7940 

0.7232 

0.6636 

0.74 

0.7102 

0.7865 

0.7136 

0.6765 

0.76 

0.7254 

0.7788 

0.7038 

0.6896 

0.78 

0.7409 

0.7709 

0.6937 

0.7029 

0.80 

0.7568 

0.7626 

0.6834 

0.7104 

0.82 

0.7733 

0.7541 

0.6727 

0.7303 

0.84 

0.7903 

0.7454 

0.6617 

0/7447 

0.86 

0.8079 

0.7360 

0.6501 

0.7600 

0.88 

0.8264 

0.7263 

0.6381 

0.7746 

0.90 

0.8459 

0.7160 

0.6253 

0.7906 

0.98 

0.8669 

0.7049 

0.6118 

0.8076 

0.94 

0.8897 

0.6927 

0.5970 

0.8260 

0.96 

0.9154 

0.6791 

0.5806 

0.8464 

0.98 

0,9463 

0.6627 

0.5609 

0.8706 

1.00 

1.0000 

0.6342 

0.5274 

0.9119 
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TABLE 4. q* 2 as a function of M 


M 

Ci *) 1 8 

M 

(q*) 2 

0.00 

0.0000 

0.50 

0.4391 

o.os 

0.0008 

0.52 

0.4682 

0.04 

0.0033 

0.54 

0.4975 

0.06 

0.0075 

0.56 

0.5268 

0.08 

0,0140 

0.58 

0.5561 

0.10 

0.0207 

0.60 

0.5853 

0.12 

0.0298 

0.62 

0.6142 | 

0.14 

0.0404 . 

0.64 

0.6439 

0.16 

0.0526 

0.66 

0.6712 

0.18 

0.0663 

0.68 

0.6991 

0.20 

0.0814 

0.70 

0.7262 

0.22 

0.0979 

0.72 

0.7528 

0.24 

0.1158 

0.74 

0.7787 

0.26 

0.1350 

0.76 

0.8039 

0.28 

0.1553 

0.78 

0.8280 

0.30 

0.1776 

0.80 

0.8512 

0.32 

0.1995 

0.82 

0.8732 

0.34 

0.2233 

0.84 

0.8941 

0.36 

0.2477 

0.86 

0.9136 

0.38 

0.2732 

0.88 

0.9317 

0.40 

0.2994 

0.90 

0.9483 

0.42 

0.3263 

0.92 

0.9632 

0.44 

0.3538 

0.94 

0.9763 

0.46 

0.3819 

0.96 

0.9872 

0.48 

0.4103 

0.98 

0.9954 



1.00 

1.0000 
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Figure 2 
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